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Abstract 

We study the off-equilibrium two-point critical response and correlation func- 
tions for the relaxational dynamics with a coupling to a conserved density 
(Model C) of the 0(N) vector model. They are determined in an e = 4 — d 
expansion for vanishing momentum. We briefly discuss their scaling behav- 
iors and the associated scaling forms are determined up to first order in e. 
The corresponding fluctuation-dissipation ratio has a non trivial large time 
limit in the aging regime and, up to one-loop order, it is the same as that of 
the Model A for the physically relevant case N = 1. The comparison with 
predictions of local scale invar iance is also discussed. 
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I. INTRODUCTION 



Non-equilibrium dynamics of statistical systems is currently under intensive theoretical 
investigation, and new dynamical behaviors have been recently discovered in models of 
disordered systems. One of the most striking of them is aging, i. e. a persistence of the 
system in a non-equilibrium state even after a macroscopic time has elapsed since the latest 
perturbation acting on it. As a consequence, there is no "memory loss" of the thermal 
history of the system and its response to an external field, for example, will depend on it. 
This fact is commonly observed in glassy systems [1,2]. It has been pointed out [3], however, 
that this kind of behavior may be also observed in critical non-disordered models. In these 
cases the presence of slow-relaxing modes could keep the system in a non-equilibrium state 
even asymptotically for large times. Consider, indeed, a system in a generic configuration 
and, at time t — 0, bring it in contact with a thermal bath at a given temperature T. The 
resulting relaxation process is characterized by a transient behavior with off-equilibrium 
evolution, for t < tr, and a stationary equilibrium evolution for t > Tr, where Tr is the 
relaxation time. In the former the behavior of the system is expected to depend on initial 
conditions, while in the latter time homogeneity and time reversal symmetry (at least in the 
absence of external fields) are recovered and such a dependence is lost; fluctuations are thus 
described in terms of "equilibrium" dynamics. 

In the following we focus on ferromagnetic systems </?, quenched at their critical tempera- 
ture T c [4] for t — (interesting behaviors are observed also in the case of non-instantaneous 
quench, i. e. for time dependent thermal bath [5]). A convenient way of describing dynamics 
is to study two-time response and correlation functions. The former is usually defined as 
Rx(t, s ) = ^(v 9 x(^))/^^(s) ) where f is the magnetic order parameter, h is a small external 
field applied at time s > in the point x = 0, and (•) stands for the mean over the stochas- 
tic dynamics. The latter, instead, is defined as the order parameter correlation function 
C x (t,s) = (<p x (t)<p (s)). 

If the system does not reach the equilibrium, the response and correlation functions will 
depend both on s (the "age" of the system, also called "waiting time" ) and on the observation 
time t. To characterize the distance from equilibrium of an aging system, evolving at a fixed 
temperature T, the fluctuation-dissipation ratio (FDR) is usually introduced [6,3]: 

TRJt,s) 

X * {t - s)= £m^>- (L1) 

When t and the waiting time s are both greater than tr, the dynamics is homogeneous in 
time and time-reversal invariant so that the fluctuation-dissipation theorem can be applied, 
leading to X x (t, s) — 1. This is no longer true in the aging regime [3]. It has been argued 
that the long time limit of the FDR at criticality 

X°° = lim limX x=0 (t,s) , (1.2) 

is a novel universal quantity of non-equilibrium critical dynamics [7-9]. Correlation and 
response functions have been exactly computed for a Random Walk, a free Gaussian field, 
and a two-dimensional XY model at zero temperature and the value X°° = 1/2 has been 
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found [3]. In the case of the d- dimensional Spherical Model [8], one dimensional Ising- 
Glauber chain [10,7] and two- and three- dimensional Ising Model, investigated by Monte 
Carlo simulations [8], X°° has values ranging between and \. 

Field-theoretical methods have been proven a powerful tool for the computation of uni- 
versal quantities (such as critical exponents) in critical phenomena (for an updated review 
see Ref. [11]). In this framework the problem of critical relaxation from a macroscopically 
prepared initial state has been analyzed since some years, and a new universal exponent 
associated with it has been introduced as a consequence of an additional time-surface renor- 
malization [12]. 

We would take advantage of these previous works to compute the critical FDR and the 
associated universal scaling functions for mesoscopic models of dynamics, overcoming most of 
the analytical difficulties encountered in the exact solutions of models with aging dynamics. 
In Refs. [13] and [14] this problem has been addressed for the dissipative dynamics (Model 
A of Ref. [15]) of the O(N) ferromagnetic model, whereas the purely dissipative dynamics 
of the diluted Ising model has been analyzed In Ref. [16]. Here we consider the O(N) model 
dynamically coupled to a conserved density (Model C of Ref. [15]). Physical realizations of 
this models are, e.g., intermetallic alloys [17], adsorbed layers on solid substrates [18] and 
supercooled liquid [19]. Also the deterministic microcanonical (p 4 model [20,21] is believed 
to be in the Model C universality class since the order parameter is coupled to the conserved 
energy [22]. 

The paper is organized as follows. In Section II Model C is introduced and its scaling 
forms are discussed. In Section III we derive the first order contribution in an e-expansion 
to the response and correlation functions for all values of s and t and we derive the FDR up 
to the same order. Finally in Section IV we discuss our results stressing their relevance for 
the issue (of applicability) of local scale invariance. 



II. MODEL C 



Let us consider the relaxational dynamics of an iV-component field y?(x, t) coupled to 
a noncritical conserved density e(x, t). This system may be described by means of the 
following coupled stochastic Langevin equations (Model C of Ref. [15]) 



^(x,t) = -fi^l| + e(x,t), 
d(p{-x,t) 



Se(x, t) 



+ CCM) , 



(2.1) 
(2.2) 



where 7i[ip,e] is the Landau-Ginzburg Hamiltonian for the fields (p and e with a coupling 
term between them 



H[<p,e] = J d d x 



1 /-— . \ 9 1 2 1 4- ^ 2 ^ 2 

^(V^) +2 r oV? +^9of 



(2.3) 



where Q and p are the kinetic coefficients, r oc T — T c , g and 7 the bare coupling constants, 
£(x, t) and C( x ^) zero-mean stochastic Gaussian noises with 
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(&(x,0&(x',O) = 21U(x " *)S(t - t')S ij: 

(C(x, t)c(x', o> = -2p n v 2 5(x - *)6(t - 1') 



(2.4) 
(2.5) 



The coupling between e(x,t) and <p(x.,t) does not change the static properties of the latter 
as it can be seen by computing the effective Hamiltonian for the ip field (see Ref. [23]). 
Moreover e-field static correlation functions are related to p 2 -G.eld correlation functions. 

Dynamical correlation functions, generated by the Langevin equations (2.1) and averaged 
over the noises £ and (, may be obtained by means of the field-theoretical action [24,23] 



S[(p,(p,e,e] = J dt J d d x 



(fO t (f + ll(f- 



5ip 



6e 



(2.6) 



where <p(x.,t) and e(x, t) are the response fields associated with </?(x, t) and e(x, t), re- 
spectively. It is easy to read from Eq. (2.6) and (2.3) the interaction vertices, given by 
—Qgo<pip 3 /3\, as in the case of Model A, — n^yeipip and pQ'-f <p 2 V 2 e/2. 

In Ref. [12,25] this formalism was extended to deal with relaxation of the system from 
a macroscopically prepared initial state. To take into account the effect of such initial 
condition on the dynamics described by Eq. (2.6) one has also to average over the possible 
initial configurations of both the order parameter v 9 o( x ) — v( x >^ = 0) an d the conserved 
density £o( x ) — £ ( x , t = 0) with a probability distribution e~ Ho ^ ' £o ^ given by [25] 



H [<Po\ 



d d x 



|(v. (x)-Mx)) 2 + ^-(e (x)-t;(x)) 5 



(2.7) 



This specifies an initial state u(x) for (p(x.,t) and t>(x) for e(x, t) with correlations propor- 
tional to Tq -1 and c , respectively. Response and correlation functions may be obtained, 
following standard methods [24,23], by a perturbative expansion of the functional weight 
e -{S[<p,<p,e,e\+H [(po,eo\) _ j n iti a l condition with long-range correlations may lead to a different 
universality class, as e.g. shown for the ci-dimensional spherical model with non- conservative 
dynamics [26]. 

The propagators (Gaussian two point correlation and response functions) of the resulting 
theory are [25] 



(<&(q, s)pj(-q, t)) = SijR (t, s) = 8ij 9{t - s)G(t - s) 



(ipi(q, s)ipj(-q, t)) = 6ijC9(t, s) 



Sij 



G(\t-s\) + 



r + q' 



-l)G(t + s) 



where 



and 



G(t) = e - n{q2+r ^\ 



(e(q, s)e(-q, t)) = R% q {t, s) = 6(t - s)G £ (t - s), 

(e(q, s)e(-q, t)) = C£,(i, s) = G E (\t - s\) + (c - l)G E (t + s), 



(2.8) 
(2.9) 



(2.10) 



(2.11) 
(2.12) 
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with 

G e (t) = e - pn{q2+ro)t . (2.13) 

As in the case of Model A and Model B, it has been shown that Tq 1 is irrelevant (in the 
renormalization group sense) so that we set tq 1 = [12,25]. 



A. Scaling forms 

When a ferromagnetic system is quenched from a disordered initial state to its critical 
point, the correlation length grows as t 1 ^, where z is the dynamical critical exponents [15] 
and t the time elapsed since the quench. So in momentum space, applying standard scaling 
arguments, the universal two time (s, t) response and correlation functions depend only on 
the two products q z t and q z s, where q is the external momentum. 

In particular general renormalization group argument suggest the scaling forms [12,25] 

OR q= o(t, s) = A R (t - s) a (t/s) e F R (s/t) , (2.14) 
C q=0 (t, s) = A c s(t - s) a (t/s) e F c (s/t) , (2.15) 

where -R q (t, s) and C q (t, s) are the Fourier transforms (with respect to x) of i? x (t, s) and 
C x (t, s) respectively, a = (2 — 77 — z)/z [27], and 9 is the initial-slip exponent of response 
function [12,25]. The functions Fc(v) and Fr(v) are universal provided one fixes the nonuni- 
versal normalization constant A R and Ac to have i^(0) = 1. 

In [13] the following quantity, related to the FDR, was introduced in momentum space 

It has been argued that the zero-momentum limit 

X™o= lim hm X«=o(t,s), (2.17) 

1 s — >oo t- — *oo 

is equal to the same limit of the FDR (1.2) for x = 0, i.e. X%L Q = X°° to all orders [13]. This 
fact allows an easier perturbative computation (in momentum space) of the new universal 
quantity X°°. Combining scaling forms and previous definitions, we find 

In recent works the notion of local scale invariance has been introduced as an extension 
of anisotropic or dynamical scaling (see Ref. [28] and references therein). Assuming the 
covariance of the response function under a suitable subgroup of the constructed group of 
local scale transformations, it has been argued that [28] 

R*(t, s) = i? x=0 (t, s)$(|x|/(* - s) 1 ^), (2.19) 

where [29] 

Rx =0 (t,s) = A R (t-s) a '(t/s) e , (2.20) 

and $(w) is a function whose explicit and convergent series expansion is known [28]. Fourier 
transforming Eq. (2.19) and setting q = one could obtain the strong prediction F R (s/t) = 
1. For the correlation function and its derivative no analogous result exists. 
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III. ONE-LOOP FDR 



In this Section we compute the non-equilibrium critical two-point response and correla- 
tion functions for the Model C up to one-loop order, for vanishing external momentum. We 
use here the method of renormalized field theory in the minimal subtraction scheme. The 
breaking of time homogeneity gives rise to some technical problems in the renormalization 
procedure in terms of one-particle irreducible correlation functions [12] so our computation 
is done in terms of connected functions. 

At one-loop order we have to evaluate, taking also into account causality [24], the ten 
Feynman diagrams depicted in Figure 1, three for the response function ((Ri), (R2) and 
(R 3 )) and seven for the correlation one ((Ci a;b ), (C 2a ,&), (C 3a ,b), and (C 4 )). 

In terms of them we have 

R g (t,s) = R° q (t,s) - ^t^Ri) + fi 2 7 2 (R 2 ) +P^V(R 3 ) + O(s 2 ,so7 2 ,7 4 ) , 

C q (t, s) = C% s) - ^±4 [( Cla ) + (C 16 )] + fi 2 7 2 [(C 2a ) + (C»)] (3.1) 

+ p^V [(C 3a ) + (C 36 )] + pfiV(C 4 ) + O(g 2 , g ol \ 7 4 ) • (3.2) 

In order to evaluate the FDR at criticality we set, in this perturbative expansion, r = 
(massless theory). We also set Tq -1 = 0, since it is an irrelevant variable [12], and Q = 1 
to lighten the notations. The first step in the calculation of the diagrams is the evaluation 
of the critical "bubbles" B lc (t), B 2c (t',t"), B 3c (t',t"), and B 4c (t',t"), i.e. the one-particle 
irreducible parts common to diagrams depicted on the first, second, third and fourth line of 
Figure 1, respectively. We have, in generic dimension d [13] 

Blc{t) -J^Y qi , ) ~~d/2-l WW 1 " * ' (3 ' 3) 

where N d = 2/(47r) d / 2 r(d/2). Given we are interested in the value of the FDR (2.16) for 
q = we evaluate, in the following, all diagrams for vanishing external momentum. Then 

for B 2c (t, s), B 3c (t, s) and B 4c (t, s) we have 

d d q 



B2 C (t,s) = j—^- d R q (t,s)Cl q (t,s) 



(2tt 

0(1 - s)[47rft(l +p)Y d/2 (/ - -s) Oo - I )(/ hsY 1 . (3. I ) 

d d q 



= 6(t-s)[Antt(l+ P )Y d l I [(t-s)- d l I -(t + Ks)- d l I \ , (3.5) 
B, c {t >8 , 8 ) = JJ0- i Cl(t, s) C% s) = ^T(d/2 - l)[fl(l + p)] 1 ^ 2 • 

. {(i_ s )l-^_ ( i + Ks) l^/2 + (co _ l)[(t _ Ks) l- rf /2_ (t + s) l- d /2| ) (36) 

where k = (1— p)/(l+p) < 1 (given that, for Model C to make sense, p > 0). Expression (3.6) 
for B 4c (t, s) is valid only for t > s, that for s > t is easily found, given the symmetry property 
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t' 

(Ri) 




s t' \ It" t 



(R 2 ) 




(c 2a ) 





(R 3 ) 



-> . 




s t' \ t 

(c 3o ) 




(C 36 ) 




(C4) 

FIG. 1. Feynman diagrams contributing to the one-loop order-parameter response ((Ri), (R2), 
(R3)) and correlation function ((Ci aj b), (C2 ,b)j (C3 0) &), (C4)). Response functions are drawn as 
lines with arrows going from the early time to the later one, whereas correlators bear no arrow 
lines. Solid (dotted) lines refer to the order parameter ip (to the conserved density e). 



B4c(t, s) = B4c(s, t). Once critical bubbles have been determined, it is easy to compute each 
diagram in Figure 1. 

Performing the required integrations and expanding in powers of e we find, for the bare 
response function, 



R q=0 (t,s) 



2f » 2 -i + 1 + 



1 + pe 



_ N + 2 _ 2 1 + p 2 - c 
2p(l + p) 



go 



24 



7o 



. t 
In - + 

s 



7o 2 



1 + P 



]n[n(t - s)] - 70 c 



1 



In 



1 — nv 



, +i 2 K(s/t;p) + 

1 — p 2 1 — re 

1 r\l 2-2 ~ ~4 ~2~ ~2\ 



where 



TZ(v;p) 



1-p 2 

and for the correlation function 



p 1 + kv 1 1 — KV 
In h In ■ 



1-p 2 2p 1 + p 



CF0(M) = _^£i + 2(is ( 1 + ^ + 



1 + p e 



12 



_ iV + 2 



0o- 



24 



-7o 



; 1 + P 2 - Cp 

2p(l+p) 



(3.7) 



(3.8) 



ln-+ 

s 



7 



-7 2 



70 ln[fit]+7o 2 



P 



- C0 C4; P ) - c C 2 (p) - ^±4 + C 2 (p) + -C 2 (-) 



+C 1 ( S 1 ;p)-C 1 (-- t ;p) 



^ + O(e 2 ,gleg ,f,fg ,ef) , (3.9) 



where we assumed t > s and we introduced g = N d g , 7 = N^jo and the functions 

c '<" ; "> = 2^7) ln|1 + - ~ KV] • (310) 

The first one is defined for— 1 < v < 1/k and p ^ 1 (we are interested only in the case 
p > 0). We note that contributions to Eq. (3.9) coming from C« are regular in the limit 

P ■ 1- 

The previous expressions for R q= o and C g =o have simple poles in e, so renormalization of 
the bare parameters is required. We use the minimal subtraction scheme in order to render 
renormalized quantities finite for e — > 0. At one- loop order it is sufficient to perform the 
following renormalizations [30,23] 

laZPZ ,with^l-^i + 0(^W), (3,2) 

to render two-point functions finite, since Z — 1 + 0(g 2 ) as known from statics [23]. 

Let us briefly recall the scenario of fixed points for out-of-equilibrium Model C [30,23,25]. 
The fixed-point values for the couplings g and 7 are determined only by the statics. We 
have g* — g* A + 67 2 *, where g* A = 6e/ (N + 8) + 0(e 2 ) is the fixed-point value of the coupling 
constant for Model A [23]. 

The value of 7 at the infrared stable fixed point depends on the sign of the specific-heat 
exponent a: 

~ 2 * _ f » stable for a < , case (I) , 

7 " \ ]v^f8) e + °( e2 ) > stable for « > , case (II) , [6 6) 

in the case (I), the dynamics of the conserved density decouples from that of the order 
parameter and we get back to Model A (at least asymptotically). At the leading order in 
e-expansion we have, for the O(N) model [11], 

4 — iV 

Q '=2(FTT + 0(f) ' (3 ' 14) 

thus the truly Model C dynamical fixed point is stable for iV < 4+0(e). In three dimensions, 
numerical calculations shows [11] that a is negative already for iV = 2, so that the Model 
C dynamics may be realized only for the three-dimensional Ising model (N = 1) that has 
positive a [11] (the two-dimensional Ising model has a = and the values of dynamical 
critical exponents for Model C are still debated [31]). 

As far as p is concerned we have three possible stable fixed points determined by equi- 
librium dynamics [30] 
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(a) p* = oo, stable for N > N^e) = 4 - [15/4 + 3/2 log(4/3)]e + 0(e 2 ); 

(b) p* = 2/N - 1 + 0(e), stable for iV < 2 + Ce\ ln(e)| and for N 2 < N < N u where 
AT 2 ( e ) = 4 - [7/2 + 3 ln(4/3)]e + 0(e 2 ); 

(c) p* = 0, which governs the critical behavior in the complement of the two regions, but 
it is a peculiar limit [23]. 

Finally, regarding the out-of-equilibrium dynamics, it has been shown that, whenever a > 0, 
the fixed point value for c is c* = [25]. 

We focus our attention on the only relevant stable fixed point of the model, i.e. (lib), 
for which 

°'=NW+8) e + 0{e2) - <3 ' 18) 

Taking into account scaling forms (2.14) and (2.15), we find the well-known critical 
exponents for Model C [23,25] (some of these results have been corrected at two-loop order 
in Ref. [32]) 

° ~ g ^2f~ 1 2p*(l + p*) +0(0 - 2(N-2 ){ N + sf + 0(e), 

and the scaling functions Fr and Fc are easily identified in Eqs. (3.7) and (3.9) with c* = 
F R (v) = 1 + r 2 [n(v, P*) - n(0; p*)\ + 0(e 2 ) , (3.17) 



F c (v) = 1+ 7 



~2* 



' hi(l-v)+C 1 (v;p*)-C 1 (-v;p*) 



l + p 



+ 0(e 2 ). (3.18) 



In particular substituting fixed point values, we obtain the scaling form (we remember that 
< v < 1 and iV < 2 + 0(e) so that no worries about the sign in the argument of the second 
logarithm arise) 

4 — N 

Fr{v) = 1 + 4(7V + 8)(iV-l) e [{N ~ 2) ln(1 + {N ~ 1)V) + NH1 ~ {N ~ 1)V)] + ° {e ) ' 

(3.19) 

and, for the physically relevant case of iV = 1, 

F R (v) = 1 - e V - + 0(e 2 ) , (3.20) 

that displays a correction to the mean-field value already at one-loop order (at variance with 
Model A [13]). 

We are now in the position to evaluate the FDR for Model C. We first note that its 
Gaussian expression is the same as that of Model A as far as (p and (p are concerned and of 
Model B (with some straightforward changes due to noncritical behavior of the conserved 
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field) for e and e. In order to evaluate the <^-field FDR we compute the derivative with 
respect to s of the two-time correlation function and consider its ratio with the response 
one: 



X q =o{t,s) = l + g — A I - 7 



2" q 



24 



1 



P 



*2 



_ i n + c 2 (p*) + -C 2 (-) + — ^ 

p* 1 + p* p* p* 2p*(l+p*) 



+C 1 (^;p*)-C 1 (-^;p*)-7^,p*) + ^ 



+ 0(e 2 ), (3.21) 



where <9i stands for the derivative with respect to the first argument. Note that, at variance 
with the one-loop FDR of Model A, this result depends on s/t. 

In the limit t — > oo, s fixed, we find an s-independent expression, 



l+g* 



N 



24 



2p* 



;i + p*)(i -p*) 2 



In 



1+p*) 2 1+p* 



4p* 



2p* 



+ 0(e 2 ) . (3.22) 



Taking into account the fixed point values of couplings we find 



1 

2 



1 + 



4-N 
N(N + 8)' 



(A-N)(2-N) A(N-iy [ { )l 



+ 0(e 2 ) 



(3.23) 



For N = 1, which is the physically relevant case into which Model C is non trivial, the result 
is exactly the same as in Model A, Af = = 1/2(1 — e/12) + 0(e 2 ), i.e. the presence of a 
coupled conserved density does not affect the value of X^ , at least up to one-loop order. 

In the e-expansion for N > iVi(e) = 4 — [15/4 + 3/2 log(4/3)]e + 0(e 2 ) a fixed point with 
p* = oo governs the critical behavior of the systems, but it probably disappears in three 
dimensions. In this limit we find (considering always a > to ensure Co = 0) 



1 + ~9*a : 



N 



N 



24 



Y* + 0(2-loop) 



(3.24) 



Once again, as it happens in all the models that have been considered so far in the literature, 
the loop corrections lead to an FDR that is less than the mean field value 1/2. 



IV. CONCLUSIONS 

In this work we considered the off-equilibrium properties of the iV-vector model coupled 
to a conserved energy density (Model C) in the framework of the field theoretical e-expansion. 
We computed up to the first order in e the critical FDR as a function of the waiting time 
s and of the observation time t. In the long-time limit, for the physically relevant case of 
one component (Ising model) this ratio has the same value as in purely dissipative Model 
A. Higher-loop calculations may clarify whether this is only a coincidence at one-loop or it 
is a deeper property. 

We also obtain the 0(e) expression for the response and correlation function for vanish- 
ing external momentum. In both the cases we found corrections to the mean-field forms. 
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Thus the result for the response function apparently disagrees with the prediction of local 
scale invariance (see Ref. [28] for an exhaustive introduction), i.e. F R (v) = 1. The same 
disagreement was already noted at two- loop order in the response function of Model A [14]. 
In that case, however, the presence of a very small prefactor in the correction makes very 
hard the detection of this effect both in experiments and in Monte Carlo simulations. In the 
present case, instead, for TV = 1, the correction in 

F R (v) = l-^ + 0(e 2 ), (4.1) 

should be large enough to be detectable. A Monte Carlo simulation of Model C could 
be helpful to clarify the nature of this disagreement. In Ref. [14] we stressed that some 
problems in the comparison between Eq. (4.1) and the predictions of local scale invariance 
Eq. (2.19) could be connected with the Fourier transformability of R x (t,s), concluding 
that some insight could be obtained by looking at the full q dependence of R q (t,s). This 
dependence was too cumbersome to be carried out in the two-loop computation of Model A 
and it is still a difficult task for the Model C dynamics at one-loop. In a forthcoming work 
we analyze the full q dependence of R q (t, s) for a ip 3 theory, showing that no problem arises 
with the Fourier transform. The nature of such a disagreement should be probably found 
in the limits of applicability of local scale invariance. 
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